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1 Introduction 



Symmetry principle is believed to be one of the important guiding principles in constructing the 
new physics models. As is well known in the context of standard model of particle physics, 
some non-trivial dynamics could be unified in the form of non-Abelian symmetry, such as 
gauge symmetry for electroweak and strong interactions or the approximate global flavor 
symmetry of quarks. In contrast, the symmetry-unified dynamics in the condensed matter 
physics is less explored and appreciated. Despite that, there was an exceptional S0(5) model 
proposed by Zhang [1] as the unified model of d-wave superconductivity (d-SC) and anti- 
ferromagnetism (AF). However, in order to explain the phase diagram of high temperature 
(high Tc) superconductivity one needs to add the explicit symmetry breaking terms in this 
model. 

In the context of AdS/CFT correspondence, the global symmetry of the boundary CFT is 
dual to a gauge symmetry of the bulk gravity. For example, in the study of meson dynamics 
of the holographic QCD, the flavor symmetry of QCD is dual to the gauge symmetry on the 
probe mesonic branes [2|. Similarly, in the recent proposal of holographic superconductor 
[3l H], the global U(l) symmetry for the charge conservation is again dual to a U(l) gauge 
symmetry in the bulk gravity. These gravity models indeed capture the essential feature of 
the dual CFTs. 

The gauge symmetry usually constrains on the dynamics more than the global one can do. 
In the context of AdS/CFT, this could imply that the constraint on the dynamics due to a 
speculated bulk gauge symmetry might uncover some emergent IR phenomenon of the dual 
CFT. One can then imagine that some originally disjoint flavor symmetries could be unified 
into a non-Abelian one at low energy due to the direct or indirect couplings among the various 
flavor currents. A typical example is the aforementioned SO (5) superconductor, in which the 
symmetry could be thought as an enlarged emerging symmetry at low energy by merging AF's 
S0(3) and d-SC's U(l). Motivated by this model and the emerging symmetry principle, it is 
then tempting to formulate a holographic model of high Tc superconductivity based on some 
underlying non-Abelian gauge symmetry. 

Naively, we can consider the gauge fields and the fundamental scalar under SO (5) gauge 
symmetry in the AdS-Schwarzschild background as the holographic dual to Zhang's unified 
theory of high temperature superconductors. We may then ask if it is possible for the black 
hole to grow the non-Abelian hairs by tuning the asymptotic values of the gauge fields, i.e., the 
corresponding chemical potentials for the dual flavor charge carriers. If so, one may wonder if it 
is possible to reproduce the peculiar phase diagram of the high Tc superconductors. However, 
a quick thought will turn down the proposal. This is because the phase diagram of the high Tc 
superconductor shows the competition between the anti-ferromagnetic order and the d-wave 
superconducting order, which is in conflict with the picture of coherent orders dictated by the 
underlying gauge symmetry. Here, by coherent orders we mean that different order parameters 
will influence each other to condense at the same temperature. This is in clear contrast with 
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the competing order phenomena in high Tc superconductivity. Indeed, in [T] one needs to 
add the exphcit S0(5) breaking terms in order to achieve the phase diagram with competing 
orders. Translated into the gravity dual picture, one needs to break the SO (5) gauge symmetry 
explicitly, which will usually lead to inconsistence as for massive gauge theories not via Higgs 
mechanism. 

We will not consider the complicated broken S0(5) case. As a first step, we beheve that 
the superconductivity with coherent orders is also an interesting physical phenomena to look 
into. We can introduce a non-Abelian gauge symmetry in the bulk to describe coherent 
orders of the boundary field theory. The coherent orders are arising from the condensation 
of the different kind of charge carriers. In the holographic QCD, these non-Abelian gauge 
fields are well known to be holographically dual to the quark or meson flavors. Similarly, we 
will interpret different component of the non-Abelian gauge field to be the holographic dual to 
some current associated with the different band carriers. From the usual AdS / CFT dictionary, 
asymptotic boundary values of the gauge fields will be dual to the chemical potentials of the 
corresponding band carriers. Unified non-Abelian symmetry can be understood as an emergent 
global symmetry due to the nontrivial interactions among the different band carriers. This is 
the origin of the coherent orders. 

It is interesting to mention that there indeed exist the real multi-band superconductors 
such as Magnesium diboride {MgB2) [5l|6] and the recently discovered iron pnictides [7]. The 
MgB2 does show the coherent orders in its superconducting ground state. More specifically, 
it has two energy gaps for two different band carriers with the same critical temperature Tc. 
However, there is no clear physical understanding why the coherent orders occur. Our study 
provides a clue to explain it by the underlying non-Abelian symmetry. 

In this paper, we will study the most simplest non-Abelian symmetric holographic multi- 
band superconductor, namely the model based on a bulk SO (3) gauge symmetry. The appear- 
ance of the multiple energy gaps with the same is an interesting outcome of this paper. Our 
results also show that a sub-sector of this model reproduce the coherent orders of the 2-band 
superconductor. This may imply that the underlying dynamics of MgB2 superconductor has 
a hidden SO (3) symmetry at low energy. Beside these, the phase diagram for the full 3-band 
case also shows interesting feature. At this point we would like to remind the readers that our 
model is different from holographic p-wave superconductors considered in [8], where only the 
non-Abelian gauge fields are introduced, not the fundamental scalars. 

We would like to emphasize that the motivation for the holographic study of the con- 
densed matter systems is to uncover the universal behaviors of the systems, especially the 
ones dictated by the underlying symmetries. For examples, the U(l) symmetry is important 
for the holographic s-wave superconductor pj to have the mean-field like behavior, the SU(2) 
symmetry is crucial for the holographic p-wave superconductor [S], and the near horizon 1-1-1 
conformal symmetry can explain the peculiar behavior of the holographic non-Fermi metals 
P]. Again our holographic S0(3) model can explain the appearance of the coherent orders. 
This is also beyond the scope of the holographic U(l) model. Therefore, our result implies that 
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the appearance of the coherent orders are universal in the system with non-Abehan symmetry 
among different band carriers. 

The coherent orders in the real multi-band materials such as MgB2 are believed to be due 
to the phonon coupling between the different band carriers, and can be obtained from the first 
principle calculation based on some microscopic models [10] . However, in this kind of approach 
it usually needs fine-tuning of the model parameters, and lacks the physical understanding of 
the underlying dynamics. In this regard, our SO (3) holographic model provides the physical 
insight for the underlying dynamics of the coherent orders, namely, it is dictated by the 
emergent SO (3) symmetry among the different band carriers. Though our result holds for the 
strong dynamics based on its gravity dual, it should also hold for the weak coupling regime 
as long as the same SO (3) symmetry emerges there. The strength of the coupling will only 
affect the mean-field behavior quantitatively, but not qualitatively. Indeed, we will see this is 
the case as in the holographic U(l) superconductor. Our finding may trigger in searching for 
the new strongly correlated materials which are unconventional superconductor with multiple 
order parameters and have the same critical temperature. 

Another interesting challenge for our proposal is how to identify the physical electromag- 
netic U{1)em in order to calculate the holographic conductivity. Naively, the U{1)em could 
be the Cartan sub-algebra of SO (3). However, there is no unique choice since the arbitrary 
proper linear combination of the SO (3) generators will play the equivalent role. We need to 
find some additional criterion for such an identification. In analogy to QCD, the flavored SO (3) 
current is usually different from the electric current, and the latter should be invariant under 
the rotation associated with the flavor symmetry. In this way, we will identify the bulk field 
which is dual to the U{1)em, and then evaluate the corresponding holographic conductivity. 

The paper is organized as follows. In the next section we will pull out the equations of 
motion for the gauge fields and fundamental scalars, and give proper holographic interpreta- 
tion. In section 3 we numerically solve the equations of motion for the background fields in the 
probe limit, and display the phase diagrams for the holographic multi-band superconductor. 
In section 4 we identify a gauge-invariant Cartan sub-algebra as the physical U(l) coupled to 
photon, and evaluate the corresponding holographic conductivity. Finally we briefly conclude 
our paper in section 5. In Appendix, we give the numerical results for the SO (3) conductivity 
matrix. 

2 SO (3) in AdS-Schwarzschild background 

As we have already mentioned in the Introduction, we will consider SO (3) gauge fields and 
fundamental scalars in the AdS-Schwarzschild black hole background. The scalars are the 
holographic duals to the superconducting order parameters of the boundary theory. The 
dual boundary global symmetry can be thought of as enlarged unified symmetry of multiple 
U(l) order parameters of some superconductors with multi-band carriers, e.g., the 2- or 3- 
band superconductors. Microscopically, the unification of the symmetry could arise from 
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the indirect interaction among the different band carriers via the phonon couphng. As we 
will see, our model reproduces the coherent feature of the order parameters for the 2-band 
superconductors, and this may justify the hidden SO (3) symmetry of the underlying unified 
dynamics for the 2-band carriers. 

The action for our holographic multi-band superconductor model is 

S = J d'x^g{R + - - - ^V') , (2.1) 

where the scalar is in the fundamental representation of S0(3), i.e., = (na, n2, ni)"^, and 
the covariant derivative D^ip := d^(f) — iqA^cf). q is the Yang-Mills gauge coupling parameter. 
The gauge connection is in the adjoint representation, i.e., 

/ -Al -Al \ 3 
A, = ^\ Al ~Al ] =5^4r\ (2.2) 

where the r*'s are hermitian S0(3) generators obeying the S0(3) Lie algebra, i.e., 

[r' T^] = if^^T^ , tr(rV^') = 25'^ . (2.3) 

Here f^^^^s are the structure constants of the S0(3) Lie-algebra, i.e., /^^^ = = f^^^ = 1, 
etc. The explicit representations of r^'s can be read from (2.2). The field strength is then 
given by 

F;, ^ d,A^, - d.Al + qP'^A^^Al (2.4) 

or in more compact form F^,^ = -F^,^r* = d^j^Ay — d^A^ — iq[Afj_ Ay]. 

The gauge field is the holographic dual of the current consisting of the carriers in the 
i-th band, and the scalar field n* is dual to the mean field order operator O^*^ in the i-ih. band. 
The S0(3) symmetry is the aforementioned unification of three U(l) bands of the carriers due 
to some microscopic dynamics such as phonon coupling. From the action (2.1) we can derive 
the equations of the motion. The equation of motion for is 

^ df,{V^D^'(f)) - iqA^'D^(f) - m^(f) = (2.5) 
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df.iy^F''"') + qP'^Ai^F'""' = iq[(j)^T'Dy - (Z}»^r>] . (2.6) 

To mimic the dual superconductor, we should put the probe gauge fields and scalar on a 
bulk black hole background with the standard AdS-Schwarzschild metric 

ds' = ~fir)de + — +P{dx^ + dy^) (2.7) 
fij) 
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where /(r) = ^ — y-. As usual, the temperature of the black hole is T = f^fiy^, which is also 
the temperature of the dual boundary theory. 

We now consider the probe background gauge fields as following 

A^dx^ := [Ui{ry + U2{ry + U3{ry]dt , (2.8) 

and the background scalar field configuration 

:= (n3(r),n2(r),ni(r))^ . (2.9) 

Note that the ordering is the reverse of the conventional one to make its compatible with the 
labeling for the gauge field. 

The equations of motion for nj(r) are 

n'l + + -X + p^n),^^ - yn, = , (2.10) 

for i,j = 1, 2, 3. In the above we have defined the bracket vectors as follows, 

(AB), = -{A,B2 + A2B^) , {AB)2 = A,B^,-A,B, , (AB)^ = A^B^ + A^B^ . (2.11) 

The equations of motion for Ili{r) are 

„ 2 , 2g2 d(Un)i 

n- +-n:-^(nn),^^ = o, (2.12) 



along with the first order gauge constraints due to our ansatz (2.8) 

7^l = (Hsn'g - U'^Us) - 2f{nsn'^ - ^'3^2) = 0; 

7^2 = (nin'2 - n;n2) - ij^n^^x - ^2^1) = , (2.13) 
7^2 = (Han; - n'gno - 2/(^3^ - ^'3^) = o; . 

In order to solve the equations of motion, we require the above gauge constraints to be consis- 
tent with the equations of motion. This will yield a consistent set of boundary conditions for 
the gauge fields and scalar fields at horizon. To see this, we consider the proper combination 
of the field equations (2.10) and (2.12). One can easily see that out of the combinations of 
these six equations, there are three independent equations as follows 

^'. + 7^ = ^ = 5 (2.14) 

where i = 1,2,3, and Cj's are the integration constants. So, in order to be consistent with 
the constraint equations (2.13), the above three integration constants should be zero. This 
will impose the boundary conditions for the gauge and scalar fields at the black hole horizon. 
One such a consistent set of boundary conditions would be to have vanishing gauge fields 
and regularity of the scalar fields at the horizon. Interestingly, this is precisely the set of 
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boundary condition which we usually consider in the holographic set up. With these above 
mentioned choice of conditions, we can solve the equation of motions without worrying about 
the constraint equations. 

At this point, it is important to note that we can reduce this 3-band model to a 2-band 
one by setting one of the following pairs to zero, (ni. Hi), (712, 112) or {n^, We can further 
reduce to the familiar U(l) model by setting two of the above pairs to zero. These relations 
imply that the multi-band models are deeply related to the U(l) case. We will see this is 
indeed the case by the similarity of the phase diagrams. 

On the other hand, we cannot reduce a 3-band configuration to a 2- or 1-band configuration 
by gauge transformation. The reason is that the vacuum is Higgsed and the gauge symmetry 
is broken so that the configurations connected by the gauge transformation will have different 
energies, and cannot be physically equivalent. This can be checked explicitly by showing that 
one cannot reduce the number of non-zero gauge fields and the scalars for a given configura- 
tion at the same time by gauge transformation. This holds even we just perform the global 
transformation at asymptotic infinity. 



3 Phase Diagrams 



In this section, we will solve the equations of motion (2.10)-(2.12) by the numerical shooting 
method, and find out the phase diagrams. Note that there are 6 functions to be solved so that 
we write a Fortran program of shooting method to perform such a task. 

As usual, we need to impose the boundary conditions to solve the equations of motion. 
Moreover, the chosen boundary conditions should be also consistent with the gauge constraints 



(2.13). After manipulating the combinations of the gauge constraints, we find that the con- 
sistent boundary conditions are pretty much the same as the ones for the holographic U(l) 
superconductor with vanishing gauge fields and regularity of the scalar fields at the black hole 
horizon, namely, at the horizon r = tq, Ilj = so that Ilidt has finite norm, and the equations 
of motion for rii implies = Sronym'^L'^ . 

Hereafter, we will choose L = 1 and = —2 so that rii is dual to a CFT operator O'-*-' 
with conformal dimension 1 or 2. This yields the following asymptotic behaviors at r = 00, 

(1) (2) 

= ^ + \ + ... , (3.1) 

Hi = /i, - ^ + ■ ■ • . (3.2) 

r 

From the above we can read off the properties of the dual CFT, i.e., the condensate of the 
operator (9*^*^ is given by 

(0«) = v^nf), a = 1,2 (3.3) 

with eabUp = 0. For simplicity, we only consider the case of a = 2 in this paper. The value 
of {02^) at zero temperature is the energy gap for the i-th band carrier to form the BCS-like 
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Cooper pairs, and the values of /ij and pi are the chemical potential and the carrier density of 
the z-th band carriers, respectively. 




Figure 1: ^ v.s. ^ , i = 1,2 — the phase digram for the holographic 2-band supercon- 

J-c J-c 

ductor, which is only a function of '■= We show the cases with = 1, 3/4, 1/2. 

Since the 2-band superconductor is well studied by the experiments, we will first focus on 
the the 2-band case of our model by setting 713 = Ha = 0, i.e., we will consider only the pairs 
(ni,ni) and (712,112). The phase diagram from our numerical result is shown in Fig. [Tj We 
see that the phase diagram in terms of the dimensionless quantities is universal, and is only 
function of pi/p2- More importantly, the phase diagram shows coherent orders, and each order 
parameter obeys the BCS-like universal scaling behavior, i.e., the carriers of 2 bands condense 
at the same Tc with the universal critical behavior as the real MgB2 2-band superconductor 
does [6]. By the numerical fitting, we find 

(O?) ^ 163Tf — il=(l - ^)^/^ 7 = 1, 2 for T ^ T, . (3.4) 

This is in contrast to the case for U(l) holographic superconductor, {P2) — 144T^^(1 — ^Y^"^- 
However, in both cases we all have the mean-field critical exponent (5 = 1/2. 

Moreover, the scaling relation between and the carriers' densities is as follows from the 
numerical fitting 



r,~0.118^p? + p2. (3.5) 

This is in analogy to the one for the U(l) case, i.e., Tc — O.llSp^''^. 

Up to now, our results of the phase diagram agree well with the BCS-like behavior for 
the 2-band superconductor, it could be the strongly coupled version of the ordinary 2-band 
superconductor. 



7 



Now we turn on all 3 bands at the same time, numerically we find that the results are more 
sensitive to the intrinsic numerical errors as expected. The phase diagrams show the similar 
mean-field feature as the 2-band case. Some of the typical phase diagrams are shown in Fig. 
|2] It is interesting to see that by tuning the chemical potentials one can collapse the 3 
bands into 2-band or 1-band cases. Moreover, there is a inversion of the vevs or the energy 
gaps as shown in Fig. and (c). 

10 I J 10 J I 10 I I 10 I 




0.0 0.2 0.4 0.6 0.8 1.0 

T/Tc (a) 




H=3.0 

n,=6.o 
H=3.0 



0.0 0.2 0.4 0.6 0.8 1.0 

T/Tc (b) 



0.0 0.2 0.4 0.6 0.8 1.0 

T/Tc (c) 



0.0 0.2 0.4 0.6 0.8 1.0 

T/Tc (d) 



Figure 2: v.s. z 



1, 2, 3 — the phase digram for the holographic 3-band supercon- 



ductor. We show the band-gap competition by tuning the chemical potentials. 

As noted, the temperature dependence of the {02^)^s near the critical point still conforms to 
the mean-field behavior, however, the chemical potential dependence is far more complicated 



than (3.4) for the 2-band case. We cannot find the complete dependence on the chemical 



potentials, but just the proportionality relation as follows. 



with 



Ci : C2 ■ C3 



c 



1,2 



for T ~ 



2/i2|/^3 - 


- 1^2 




1/^1 - 


1 + 







2/i2|/il - 


- At2 




1/^1 - 1^2 


1 + 


- 


- /^2| 



(3.6) 



(3.7) 



On the other hand, the scaling relation between the critical temperature and the carriers' 
densities has the similar form as the 2-band and the U{1) cases, namely. 



0.118 



pI + pI + pI 



(3. 



So far we have discussed about the superconducting phase diagram of a holographic multi- 
band superconductor. Behavior of the condensation for each band is universal in terms of 
temperature. As we decrease the temperature, condensation happens at the same temperature 
for all the band carries. As we have mentioned in the introduction this phenomena is quite 
similar to the real two-band superconductor MgB2. In the subsequent section we will calculate 
the superconducting transport properties like optical conductivity under small electromagnetic 
perturbation. 
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4 Holographic Conductivity 



In this section, we would like to derive the field equations for the gauge field perturbation 
on the above background. We then solve these equations to obtain the holographic real time 
Green functions of boundary currents, from which we can extract the conductivities. 
Let us turn on the gauge field perturbation of the x-component as follows, 

SA^dx^" := e-^('^*-^2^)[ai(r)r^ + aiir^ + as{r)T^]dx , (4.1) 



Plugging the perturbed fields into the field equations (2.6), and expand it up to the linear 
order, from which we can derive the field equations for the perturbed fields. The results are 
the following, 

<+J<+j2 (^'«^ + 2^ga;r'n,afc + q\n ■ Ua, - U,U ■ a))-i (^2q\an),^^^ + '^a^ = 

(4.2) 

— > — » 

where we define 11 := (111,112,113), a := (01,02,03), and the inner product such as 11 ■ a = 
11101 + 11202 + 1^303. In the above expressions, we have again used the bracket notation defined 



in (2.11). We also have noticed that there is no constraint equations among Oj's perturbations 



as long as we consider a specific form (4.1) of the electromagnetic perturbation. Therefore, 
from the holographic point of view the external electric field at the boundary does not produce 
any scalar perturbation in our background. 

Similar to the abelian case, we then solve the above equations by imposing the incoming 
wave boundary condition in the near horizon region, namely Oj = /^*'^'^^/^''"[l+Oi,i(r— ro)+- ■ ■ ]. 
Then from the asymptotic behavior of all the fields at the asymptotic boundary we found 

aj{r,k,u) = af\k,u) + ^{k,u) + ■ ■ ■ , (4.3) 

It then seems that one can evaluate the holographic conductivity by the standard holographic 
prescription for the Ohmic law [HI [121 UHl E]. However, one should be careful about the fact 
that all the perturbed fields Oj's are linearly coupled to each other in the bulk. This coupling 
among the fields implies that source of one particular perturbed field will also source the 
current corresponding to the other components. If we identify Oj as dual to the U(l) current 
for the i-th band, this will then yield 

{J,)=a,,E„ (4.4) 

with 3x3 matrix general conductivity matrix in the global SO (3) space. Noe that the 

i,j are the indices for the S0(3) internal space, not the space-time ones. Since the conductivity 
matrix cannot be measure directly, we will not discuss it further here but in the Appendix. 

In our model we have three different kinds of carriers corresponding to three components 
of a fundamental scalar field in the bulk. Each carrier has the corresponding boundary current 
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dual to . This dual current transforms like a vector under the global internal symmetry. This 
is, however, not the current that we are interested in. Now how do we define the conductivity 
under U(l) electromagnetic field? According to the standard holographic prescription, the 
above currents are actually some kind of weekly gauged fiavour currents similar to the standard 
s-wave superconductor. As we have argued, the system we are considering is a holographic 
multi-band superconductor. So, all the band carriers should couple to a single physical photon 
field through which we will calculate the optical conductivity of the system. 

Naively, we can identify one of the three generators of S0(3) as the physical U{1)em, 
namely the Cartan subalgebra, and the other two as the ladder operators. More specifically, 

[r^r±]=±r^, [r+,r-] = rl (4.5) 

However, the above identification is not unique, and any linear combination of the three 
S0(3) generators with proper normalization will play the same role. Even though there is an 
arbitrariness in defining the U(l) subgroup, any physical quantity should be independent of 
that particular choice. But this does not happen in our case. Specifically, considering above 
particular definition, we can find out the solution for 03 which coupled to a boundary current 
with the source of physical electric field, but other two components ai and 02 are sourceless or 
normalizable . Then one can calculate (T33 and identify it as a physical conductivity using the 
standard holographic prescription. This is how we indeed calculated the conductivity matrix 
as shown in the Fig|5]of the Appendix. We can apply the same procedure for the other choice 
of U(l), but those choices lead to different values of the conductivity as one can easily see 
from the plot for the conductivity matrix. This discrepancy is intuitively obvious because 
the background chemical potential for each band carrier is different. So, when we consider 
different choice of U(l) subgroup of the bulk gauge group, we are actually not considering 
same photon field at the boundary. 

Therefore, the Cartan subgroup cannot be the physical U{1)em for the electromagnetism. 
Instead, the physical U(l) electric current should be independent of the choice of the Car- 
tan subalgebra, invariant under the global SO (3) rotation and also does not depend on the 
intrinsic properties of the individual band carrier. We therefore can interpret that particular 
combination as a dual to the physical U{1)em current at the boundary. Accordingly, we can 
evaluate the holographic conductivity which is invariant under the SO (3) transformation and 
also does not depend on individual chemical potential. 

It seems strange that the physical U{1)em current is defined only for the perturbation but 
not for the background. However, in the holographic approach this is natural because the 
background gauge field is dual to the chemical potential and carriers' density, which define the 
properties of the vacuum not the dynamics. In fact, the different band carriers should have 
the corresponding chemical potentials and carriers' densities as specified by the background 
profiles. In contrast, the perturbation of the gauge field is dual to the electric source and the 
djTiamical current in the context of linear response theory. In this case, all the carriers are 
coupled to the same U{1)em source, and an unambiguous identification of physical U{1)em 
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Figure 3: The AC conductivity ag.^ of 2-band superconductor for /ii = 2.5, fi2 = 3.0, /is = 0.0. 
(Left): the real part of conductivity, and (Right): the imaginary part. 



IS necessary. 

Mapping the above consideration back to the bulk point of view, the SO (3) rotation is 
translated into the gauge transformation for the total gauge field, namely, 

= d,a' + qr'Aj^a' (4.6) 

where the total gauge field includes both the the background ansatz ( |2.8[ ) and the fluctu- 



ation one (4.1), and a*'s are the gauge functions. Since we are only interested in the S0(3) 
rotation relating different choices of Cartan U(l), this implies that the gauge function func- 
tions should be independent of the boundary coordinates. This particular class of gauge 



transformation retain the transformed total gauge field in the same clas^of the ansatz (2.8) 



and (4.1) if we also require the gauge functions are also independent of r. Then, the gauge 
functions become constant gauge parameters and the gauge transformations reduce to the 
global rotation, and the only nontrivial parts of the transformations are 

SUi = P%a'' , 6ai = f^^a^a''. (4.7) 

The task to find the dual field of physical U{1)em is then equivalent to find the appropriate 



linear combination of the gauge fields a^'s so that it is invariant under the transformation (4.7). 
it is then straightforward to see the field 11 ■ a satisfies this constraint, and we will identify it 
as the holographic dual to the physical U{1)em current. The peculiar dependence of this dual 
current on the background profile 11 may reflect the physical relevance that the current here 
is actually the vev of the current operator, which will then depend on the properties of the 
vacuum encoded in 11. 



^Namely, the background is a time-component of the gauge field as only a function of r, and the perturbation 
is the x-component of the gauge field. 
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Figure 4: The AC conductivity a^.^ of 3-band superconductor for /ii = 3.0, /i2 = 2.5, /is = 2.8. 
(Left): the real part of conductivity, and (Right): the imaginary part. 

Accordingly, we can define the holographic conductivity cr^j g., and using the asymptotic 



expansion (4.3) it can be expressed as 



n-a 



lim 

fc-5>0 



(0)^ 



(0) 



(4i 



Our numerical result for cTj^ ^ for the typical 2- and 3-band cases are shown in Fig. [3] and 
Fig. |4} respectively. We also consider the other cases, including the 2-band one considered 
in the last section, and the results are similar. We shall also mention that in our numerical 
calculation we find that ag.g. is independent of the choices of the initial conditions for a^'s at 
the black hole horizon. This is not so trivial since the perturbed fields a^'s are linearly coupled 
to each other in the bulk. This property thus supports the identification of 11 ■ a with the dual 
of U{1)em current is sensible according to the holographic prescription. 

In Fig. [3] and Fig. |4| we see that the gap appears clearly in the real part of the total 
conductivity Re[(Ttof] as expected from the Ferrell- Glover sum rule to make up the carriers for 
the infinite DC conductivity. From the tail near zero frequency, we can extract the density of 
states for the normal component of the carriers, i.e.. 



lim Re[atot{uj)] ~ Cexp (-7A/T) 



(4.9) 



where A := [ELi(C?)T/V2. 

Moreover, the imaginary part of the total conductivity has a pole at w = 0. From the 
Kramers- Kronig relation this implies the infinite DC superconductivity caused by the non- 
zero superfluidity density, Ug. From our numerics, we can extract the scaling behavior 



Us := hm uIm[atot{uj)] ^ D{Tc- T) 



as 



T^Tr. 



(4.10) 
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The numerical fitting gives C ~ 14, D ~ 24 and 7 ^ 0.97 for the 3-band case, this is the 
typical mean-field hke behavior as for the holographic s-wave superconductor [3J. For the 2- 
band case, we find the numerical fitting values of C, D and 7 are not universal but depend on 
the chemical potentials. However, we cannot find their dependence on the chemical potentials 
in the closed form. 

5 Conclusion 

In this paper, we study a holographic model which exhibits the low energy behavior of a 
multi-band superconductor. Specifically, the two-band superconductor like MgB2 has been 
studied quite extensively from the theoretical and as well experimental point of view. Most 
of the properties of this kind of superconducting material are believed to be explained by 
the standard BCS theory, however, the insightful understanding for the appearance of the 
coherent order is obscure in the first principle calculation based on BCS theory. In this report 
we tried to construct a holographic model of this kind of multi-band superconductor, and 
attribute the coherent orders to the underlying emergent SO (3) symmetry. We conjectured 
that the interactions among different band carriers are dictated by an underlying SO (3) global 
symmetry, which is then dual to the bulk SO (3) gauge dynamics. Surprisingly, our model re- 
produces the phase diagram with the desirable feature for multi-band superconductor, namely, 
the feature of the coherent orders. Moreover, we identify a gauge invariant linear combination 
of the perturbed gauge fields as the dual current coupled to the physical U{1)em photon, 
and it shows the mean field BCS-like behavior for the holographic conductivity. However, it 
deserves further study to clarify the physical subtlety of such an identification, especially its 
dependence on the background profile. Finally, it is a natural next step to see if these features 
remain intact after taking into account the back reaction of the bulk fields to the background 
geometry. 
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Appendix 

The general definition of conductivity matrix is 

aij{u) = -lim I ci^a;e*"*-'^"e(t)([J,(t,x), J,(0,0)]). (5.1) 

^ fc-s>0 J 
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This quantity may not be measured directly in the real experiments since it should be highly 
nontrivial to tune the electric field among different band carriers. Despite that, the conductiv- 
ity matrix provide some "microscopic" picture for the holographic multi-band superconductors, 
and it is interesting to find out its behavior. 



By fixing the background (2.8), it is easy to see that gauge field perturbation 
invariant under the infinitesimal gauge transformation, i.e., ^a, = D^a^ where = + gA^ 
is the covariant derivative with respect to the background gauge field A^, and a^'s are the 
gauge parameters. Therefore, by using AdS/CFT correspondence, one can define the S0(3)- 
invariant general conductivity matrix as follows 



aij{u) = - hm ^ 1^(0) . (5.2) 

fe-*o ujaj '[k, 00 



A typical plot for conductivity matrix is shown in the Fig. |5] for 2-band holographic super- 
conductor. Note that it does not have the mean-field like behavior. 
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Figure 5: Conductivity matrix of a 2-band superconductor for ^ = 0.8460, yU2 = 3.2, /is = 2.5 
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